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Abstract. Supemovae (SN) and supernova remnants (SNR) have key roles in galaxies, but their physical descriptions is still 
incomplete. Thus, it is of interest to study neutrino radiation to understand SN and SNR better. We will discuss: (1) The ~10 
MeV thermal neutrinos that arise from core collapse SN, that were observed for SN1987A, and can be seen with several 
existing or planned experiments. (2) The 10-100 TeV neutrinos expected from galactic SNR (in particular from RX J1713.7- 
3946) targets of future underwater neutrino telescopes. 
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1. INTRODUCTION 

In the economy of a galaxy, the stars with mass > 6 — 
10 Mq are known to play a prominent role. The last in- 
stants of their brief life (on ten million year scale) is a 
crucial moment. In fact, the gravitational collapse of the 
'iron' core leads to supernovae of type II, lb, Ic and even- 
tually to pulsars, neutron stars and stellar black holes. 
The gas subsequently expelled by all supernovas includ- 
ing la (that lead to the observed supernova remnants) 
is suspected to be the place where cosmic rays are ac- 
celerated. A complete theoretical understanding of these 
systems (and in particular of the supernova) is not yet 
available. However, it is known that neutrinos of ~ 10 
MeV are emitted during and just after the gravitational 
collapse, and that neutrinos above 1 TeV can be pro- 
duced in the supernova remnant at least in special cir- 
cumstances. Therefore, neutrinos can be used as diagnos- 
tic tools of what happens in these interesting and com- 
plex systems. With these considerations in mind, we will 
discuss here briefly neutrinos from supernova and from 
supernova remnants, aiming to outline the main 'what', 
'when', 'where', and 'how' questions. 

Outline and references. In the rest of this section, 
we discuss a number of general facts pertaining to su- 
pernovas and supernova remnants; the references that we 
use are d 01 . In the second section we concentrate on 
supernova neutrinos; we use 1 1, 2], and also |3]. In the 
third section we describe the calculations of |4|. In these 
4 papers, one can to find more details and a complete list 
of references. In the appendix, we show in details a sim- 
ple way to calculate neutrinos from cosmic rays, adopted 
also in |4]. 

Generalities. Before passing to neutrinos, we would 
like to discuss a few general questions. 



Guessing 'where' Let us recall that the only SN within 
the capability of existing neutrino detectors are those 
from our Galaxy (or some irregular galaxy around us). 
The best guess we can propose for next galactic core 
collapse SN is 

(L) =10 ±4.5 kpc (1) 

and it is motivated as follows: 

1) We are R = 8.5 kpc from the galactic center. 

2) The distribution of the matter that can go in super- 
novae is: p ~ re~ r ' r ° with r=distance from the center and 
ro = 3 kpc, possibly summing a 8 (r) to describe the 'bar' 

3) We calculate the distribution in function of L=distance 
from us, integrate over the galactic azimuth and get the 
result above. 

Guessing 'when' The rate of core collapse SN in the 
Milky Way is expected to be 

Rsn= 1/(30-70 years) (2) 

The most reliable method is: count SN in other galaxies, 
and correlate with galactic type. 1 A similar rate expected 
for SN la. Possibly, we missed several galactic SN due 
to dust. 2 For the future, we will have better coverage of 
galactic SN with neutrinos, IR and perhaps with gravi- 
tational waves. From the absence of neutrino bursts, one 
derives 

R SN > 1/(21 years) at la (3) 

We assume Poisson statistics (namely, exp(— TR$n) = 
1 — C.L) with T = 24 years. To obtain the last number, 



1 The Padova-Asiago database |5] includes several thousand SN. The 
Milky Way could be Sb or Sb/c, which means a factor 2 uncertainty. 

2 Indeed, with about 10 supernovae seen in last 2000 years, in order to 
have a SN each 25 years we need to admit we see just one SN each 8, 
that is about the right figure on accounting the presence of the dust and 
the possibility that a supernova explodes on the day sky 01 . 



we note that till 1986 only Baksan |0] worked with 90 % 
DAQ livetime, then we assume 100 % coverage. 

We conclude by listing the number of SN pre- 
cursors and descendants. We estimated these 
populations by assuming R^L = 1/(25 years). 



Object 


Lifetime 


Number 


Pre-SN with v 
Pulsars 
SNR 
young SNR 


20 million y 
2 million y 
100.000 y 
2000 y 


400.000 
40.000 

2 x 2000 
2x40 



By 'Pre-SN with v' we mean core collapse SN. Recall 
that core collapse SN produce neutron stars (NS) or 
stellar black holes, and that pulsars are 'active' NS. Type 
la SN make white dwarfs, instead, but are also supposed 
to produce SNR (which explains the factor 2 above). 



2. vS FROM SUPERNOVAE 

We begin recalling the astrophysics of core collapse. 
The giant stars-pre-SN-burn in sequence H, He, C and 
Si, Ne, Mg, Na etc, and form an "onion structure", with 
an inert 'iron' core. Violent stellar winds occasionally 
modify external mantle in latest stages; apparently, this 
happened for SN1987A that was a — 20 M Q blue gi- 
ant. What happens in the core? The gravitational pres- 
sure is balanced by e~ degeneracy pressure, and the core 
grows. As demonstrated by Chandrasekhar, the pressure 
of free electrons is unable to supply an equilibrium con- 
figuration when e~ become relativistic. Now, the iron 
core mass is ~ 1.4 M , the collapse begins and the se- 
quence of the events becomes uncertain (more on the ref- 
erence picture, the so-called "delayed scenario"). What 
about the energetic of the collapse? The total energy 
of the collapse is very large, and 99 % of this energy 
is carried away by neutrinos. With M„ s /Mp, = 1—2, 



Rn 



15 km(M Q /M m 



1/3 



, we estimate 



= (1 -5) x 10 53 erg 



7Rns 

The delayed scenario. This is a pictorial summary 
of the 'delayed scenario' by Wilson & Bethe 0]: 
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The energy radiated in any neutrino species 
e,e,lJ,,jl,X,Z is expected to be the same within a 
factor of two |9|: 

<S(> r *~* J (§g r ^ J S^x 

where x denotes any among /x, fi, x, f. Indeed, in this 
picture non-electronic neutrinos and antineutrinos are 
produced in a similar amount. 

Note that neutrinos are mostly emitted in cooling (80- 
90 %) and accretion (10-20 %) phases. The latter one is 
the phase crucial to understand the nature of the explo- 
sion, whereas the first phase is almost thermal. In this 
sense, the ignorance of the detailed mechanism of the 
explosion does not preclude a description of the bulk of 
neutrino radiation (to put it with a slogan, ignorance is 
self-consistent). 

Now we come to a prescription for the fluence: 

F .(E) - -£l ^_ 7 « e -(«+l)z z - — 

,( > ~ 4M? (Ei) 2Z ' (Ei) 

where (£;) is the average energy of the neutrino species 
i = e,e,x; N ensures that the total energy carried is 
Si. (If one needs to describe time dependent situations, 
Si -> Li{t) = dSt/dt, (Ej) -> (Ei(t)), a -> a(t).) The 
expectations for time integrated quantities are: 



(Eg) = 12-18MeV, 
S g = (2- 10) x 10 52 erg 



(E X )/(E S ) = 1-1.2 
S x /S g = 1/2-2 



One guesses S e = Sg (not so important); the v e average 
energy can be estimated from the emitted lepton number. 

Oscillations of supernova neutrinos. Oscillations of 
SN neutrinos are pretty simple to describe. To account 
for oscillations we need to assign just 2 functions, P ee 
and Peg : 



*F e — F e P ee 

- F° P 

•Fe+F„+F T 



-Ffl P lie + 



F° Pr 



and similarly for antineutrinos (we will consider only 
oscillations of massive vs). The relevant densities to 
calculate P ee and Pgg are p so i ~ 10 gr/cc (He) and p atm ~ 
10 3 gr/cc (C+O). 

There is a great interest on the unknown parameter 
C/ e 3, and supernovas offer a chance to learn more on 
that. However, we should keep in mind the uncertainties 
mentioned above. Let us try to discuss this point further. 
We begin with the formula F e = F° - P ee (F® - F® ) . If we 
have normal mass hierarchy: 



Pep — 



sin 0i3 



■0 9 l3 'large', > 1° 
■0.3 On 'small', <0.1 c 



sin z 0i2 

Now we can ask a precise question on U e j,\ 
Can we distinguish the two cases? 



TABLE 1. Depending on the phase of neutrino emission, it is possible to learn more on l/ e 3 studying 
of SN neutrinos. In the table, goods and bads of each specific phase are mentioned. (LSD refers to the 
observations of the Mont Blanc detector, that according to Imshennik and Ryazhskaya 1 10] could be 
interpreted assuming an intense neutronization phase before the main v burst.) 



Emission 


Good 


Bad 


Remarks 


cooling 


strong v radiation 


uncertainties, small effect 






accretion 


strong v radiation 


uncertainties ! 






neutronization 


clean signal 


weak v radiation 


Fx r 


-0 


neutroniz.++ with rotation 


clean and strong sign. 


uncertain!! 


F?' 


- (LSD?) 



Tabled provides a check-list. Note that if the mantle is 
stripped off till densities p > 10 gr/cc (e.g., with SN Ic) 
we have vacuum oscillations, 0.3 — » 0.6 (Selvi). This is 
rare, but not impossible. 

SN1987A. The first detection of SN neutrino is of 
epochal importance. These observations fit into the 'stan- 
dard' picture for neutrino emission (see next figure), but 
there are some puzzling aspects: 

1 . 1MB O andKamiokande-II find forward peaked 
distributions; e.g., (cos 0) are ~ 2 a above expectations. 

2. ~ 15 MeV and (£™ B ) - 30 MeV (±2.5 MeV) 
are quite different even correcting for efficiencies. 

3. The time sequence of events looks different; when 
combined not so bad (but abs. time is unknown). 

4. The 5 LSD events, occurred 4.5 hours before the main 
signal, cannot be accounted for. 

In figure ^ will stress the consistency of the 'standard' 
interpretation, but the space for non-standard ones is not 
small (not only due to limited statistics). A reasonable 
agreement with expectations is obtained if: 

{Ee)=E = 12 -16 MeV 

£= (2-3) x 10 53 erg 

zero or a few v e e — > v e e events in KII 



2) We have a new SN each t$n ~ 25 year, with about 
§ ~ 10 51 erg in kinetic energy, that is 



SN : 



= 1.2 x 10 42 erg/s 



Tot 



Comparing the two formulae, we see that if a SNR is 
able to convert a fraction /cr ^ 5 — 10 % of the injected 
energy into CR, we are home (the numbers quoted above 
shouldn't be taken too seriously, but this 40-years-old 
argument maintains its appeal). 

We proceed describing the system. In 2000 years, 
the SNR proceeds by ~ 10 pc. The density is ~ 0.2 
protons/cm 3 , too faint to permit a significant produc- 
tion of secondaries. In the galaxy, the most dense non- 
collapsed objects are the molecular clouds, whose den- 
sity can reach 10 4 protons/cm 3 . The 2 objects form: 
A COSMIC BEAM DUMP 



Molecular 
cloud 



the target 



SN shell 
the accelerator 



3. vS FROM SNR 

The leftover gas is the SNR. The kinetic energy is a 
few times 10 51 erg, which means that the gas is in free 
expansion with velocities of about ~ 4000 km/sec. There 
are various phases of the SNR, according to their age, 
and various shapes (shell, plerionic, or mixed). 

An argument by Ginzburg and Syrovatskii suggests 
SNR as main source of galactic cosmic rays (CR): 
1) The Milky Way irradiates CR. Take V C r = nR 2 H with 
R ~ 15 kpc, H ~ 5 kpc as the volume of confinement. 
Take Tcr = 5 x 10 7 years as CR lifetime in the Galaxy. 
We get the CR luminosity: 



That is an ideal configuration, since 



VcR-PcR n n 1n 4l , 

cr = = 0.9 x 10 erg/s 

?CR 



CR collisions 



high energy y 
high energy v M ,V e 



Thus, we expect to have y and neutrinos. 

There is some evidence that the system called RX 
J1713. 7-3946 is a 'cosmic beam dump'. Indeed, 

1) It is seen in X-rays, with many details 

2) It is in Chinese Annales, 393 A.D. Q 

3) A molecular cloud seen in CO and 21 cm H line 

4) And most interestingly, CANGAROO (since 2000, see 
fil ) and H.E.S.S. (since 2004, see 0) do see TeV y 
rays. 

The distance is L = 1 kpc, the angular size about 1°, the 
density of the cloud ^100 part/cm 3 . The source is trans- 
parent to gamma rays, so that the neutrino flux can be 
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FIGURE 1. Average energies and number of events for three detectors. The horizontal lines describe the measured values, the 
oblique lines describe the theoretical expectations. 



calculated easily and reliably as shown in the appendix. 
(However, no item above seems to be waterproof, what 
we discuss is just an appealing interpretation.) 
We assume that the cosmic ray spectrum: 



F p = KE 



-T 



T = 2-2.4 



interacts with the molecular cloud and is the main source 
of visible gamma rays, through the reaction p — ► n° — * 
y. Neutrinos originate from similar processes, e.g., the 
reaction p — > 7T+ — > ji + — > V e will yield V e . Using the flux 
F r = 1.7 • 10- n (£/TeV)- 2 2 TeV-'cm-^- 1 measured 
by H.E.S.S. in the range between 1-40 TeV we get: 



f* =7.3X10-* 

(tSv) 



-2.2 



TeVcm-s 



-12 



= 7.4 X 10 



*?. =4.7x10"" (^) 



-2.2 



TeVcm z s 



-2.2 



1 



TeVcm-s 
-2-2 1 

TeVcm 2 s 



3 flavor oscillations take the simplest form (vacuum aver- 
aged or Gribov-Pontecorvo) and can be included easily. 
This leads to an important modification of the fluxes. 

Finally we pass to the signals of neutrinos, v interac- 
tions are due to deep elastic scattering. The simplest and 
most traditional observable is induced muons, that can be 
correlated to the source by mean of an angular cut. It is 
important to recall that high energy are to some ex- 
tent absorbed from the Earth, and that when the source is 
above the horizon, it is impossible to see anything due to 
the background from atmospheric jj.. Along with oscilla- 
tions, these effect decrease the observable neutrino sig- 
nal. For an ideal detector, located in the Mediterranean, 



with Area=l km , Datataking=l year, E,^- = 50 GeV 
we find that the number of expected events is about 10 
(this was 30 if oscillations, absorption and ji -background 
were ignored, or even 40 if the slope was F = 2.0). 



4. DISCUSSION AND PERSPECTIVES 

We do not have a clear understanding of how SN ex- 
plode. Perhaps, because it is a very difficult problem, or 
there is some missing ingredient, perhaps the answer is 
not unique (a combination of various mechanisms?), ... 
the confusion could persist even after next galactic SN. 

vs from next galactic SN have an impressive potential 
to orient our understanding. The hypothesis of an "ac- 
cretion phase", implying 10-20 % of the emitted energy, 
can certainly be tested. SN1987A does not contradict the 
'delayed scenario' seriously (but does not help much ei- 
ther). There are chances to learn on oscillations and in 
general on particle physics. The possibility to use v e and 
neutral current events deserves consideration. 

vs from SNR are an uncharted territory. Recent results 
from H.E.S.S. motivated us to consider one specific SNR 
(but we expect new results and surprises with y rays). 
Sure enough, CR acceleration in SNR is not fully under- 
stood, and there are other possible sources for TeV v as- 
tronomy, however the number of events we found (about 
10/km 2 y) suggests the need of large exposures. It is im- 
portant to improve our theoretical tools to describe SNR 
y and v. 

APPENDIX: Cosmic rays and neutrinos. We report 
here on the details of the calculation of (anti) neutrino 
spectra from the decays of secondary particles produced 



in pp interaction in astrophysical beam-dumps (table 1 
of |4]). As a first approximation, we add the contribu- 
tions to v-emission from two chains: (1) p — » m — > 
(anti)v and (2) — > wt — > /i ± — > (anti)v, where m = 
{M: ± ,^ ± } ) V = {v e ,V e ,V M ,V jU }. 

If we assume scaling, and take a power-law proton spec- 
trum Fp(E p ) = KE~ r , we can express the (anti) neutrino 
spectra from (1) and (2) as: 



where the moments: 



A p 



ih-W'W TP (F \ 



(4) 



where y/ v indicate the neutrino emissivity coefficient, 
defined in 1 16], AX the column density traversed by the 
protons and X p the interaction length of CR. Similarly, 
the y-spectrum from p — ► 7C° — ► y is: 



AX 



F y (E 7 ) = — • \ff r ■ F p (E r ) , where Xjfy ■■ 



2Z 



(5) 



In equations |4] and |3 we assume a thin target. When we 
assume scaling, following Gaisser (see 1 17], formulas 4. 1 
and 4.2 in approximation 3.31) we have: 

dn v 



F V (1) =AX 



with: 
dn 



m+=7X+,K^ 



B„ 



-rm 



dE v 



D m + {En 



1 Y m E n 



-F p (E n 



(6) 



(7) 



where B m + is the branching ratio for m + — » /i^Vj 
r,„ = m^/m 2 , (subscript m = {7r,/T}), and Zp m + are the 
spectrum-weighted momenta, defined in 1 17]. We recall 
that in the scaling hypothesis the spectrum- weighted mo- 
menta depend only from the spectral index F. We can get 
the spectrum from the same formula, replacing m + 
with mT . Evaluating the integrals, we get: 



pW-** F ( F 



) E 

m+ = 7l+,K^ 



B„ 



(l-r. 



\r-i 



and so: 



,(2) 



B m+ — — (1 -r m ) 



r-i 



(8) 



(9) 



Similarly, for y/vu' we nave ( see fl^ll - formula 
(7.14), in which we replace § n with the fac- 
tor [AXZ pm -B m -Fp(E v )]/X p , and then divide by 
AXFp/A p ): 




(ii) 



for (anti)v e and (anti)v^ are given in table 7.3 from fl7il 
(the giviy) are given in table 7.2 from Q)- We can ob- 
tain y/^ from eq. ^| replacing Z pm - with Z pm + . In our 
approximation, the (anti)v e emissivity coefficients are 
given by eq.[lO] inserting respectively Z pm + for yr Ve , and 
Z pm - for yiy e and the appropriated moments. The numer- 
ical values of k v = (Yy + Wv^)/ Vy f° r several values of 
r, obtained using the spectrum-weighted momenta from 
figure 5.5 by 1 17], are reported in table 1 by |4J]. 3 
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